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We consider a system of non-interacting charged particles moving in two dimensions among fixed 
hard scatterers, and acted upon by a perpendicular magnetic field. Recollisions between charged 
particles and scatterers are unavoidable in this case. We derive from the Liouville equation for this 
system a generalized Boltzmann equation with infinitely long memory, but which still is analytically 
solvable. This kinetic equation has been earlier written down from intuitive arguments. 

I. INTRODUCTION 



> 

in 

The Lorentz model where a particle moves among and collides with fixed scatterers has provided a rich testing 
ground for kinetic theory. In particular, the Boltzmann equation is not only exactly solvable in this model, but 
the equation itself is exact in the Grad limit (to be defined below). The Stosszahlansatz used in constructing the 
Q\ ' Boltzmann equation relies implicitly on the assumption that the particle never returns to a scatterer after having 
ON ■ collided with it. The probability of such recollisions vanishes in the Grad limit. On this basis the Boltzmann equation 
for the standard Lorentz model is taken to be exact in the Grad limit. 

Recently, however, an interesting exception was discovered to this state of affairs. Let the particle have electric 
charge — e and move on a plane pierced by a perpendicular constant magnetic field B. The particle moves along 
circle arcs between collisions, and if it does not encounter a new scatterer along the arc, it will recollide with the 
initial scatterer. This destroys the above assumption and renders the Boltzmann equation invalid. Note that both 
the (two-)dimensionality of space and the presence of a perpendicular magnetic field are essential here. 

In Refs. a generalization of the conventional Boltzmann equation was proposed that takes into account consec- 
^ . utive recollisions with the same scatterer. The arguments leading to this equation were intuitive, on the same level as 
the Stosszahlansatz itself. The aim of the present paper is to provide a microscopic underpinning to the generalized 
Boltzmann equation in the form of a systematic derivation from the Liouville equation. As a background we start by 
summarizing the intuitive arguments leading to the generalized Boltzmann equation. 

The charged particle moves on a plane of (large) area A, with TV randomly placed hard disk scatterers of radius a. 
We denote by n — N/A their number density. The disks do not overlap. 

The generalized Boltzmann equation describes the evolution of the probability density /i(x, v, t) for finding the 
moving particle at time t at position x with velocity v. This non-markovian kinetic equation has the form 



D 
~D~t 



[t/T ] 

/ G (x, v, t)=na V e~ vkT ° / rfn(v • n)[0(v • n)6 n + 0(-v • n)] f G (x, S " fe v, t - kT ) , (I) 



where 



f g< „ ^ _ J /i(x,v,t) if < t <T 



f (x ' v ' t} "^ (l-e-^o )/l(x ; v ; t) if t >T , (2) 



*Permanent address: Keldysh Institute for Applied Mathematics, Russian Academy of Science, f 25047 Moscow, Russia 
^Permanent address: Department of Physics, Norwegian University of Science and Technology, N-749f Trondheim, Norway 



1 



and where v = 2|v|na is the collision frequency and To the cyclotron period. Furthermore, 

(3) 
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is the generator of free cyclotron motion with frequency |u>| = w = eB/m, and [t/T ] the number of cyclotron periods 
T = 27t/oj completed before time t. The angular integration over the unit vector n in (|l]) is over the entire unit 
sphere S 1 centered at the origin. In the gain term (positive contribution), there appears the distribution acted upon 
by the operator 6 n , defined by 

M(v) = 0(v - 2(v • n)n) , (4) 

where <fi is an arbitrary function of v. The precollisional velocity v' = v — 2(v-n)n becomes v after the elastic collision 
with the immobile (infinitely massive) scattering disk. Note that v' • n < 0. In the loss term (negative contribution), 
the precollisional velocity, v, is also from the hemisphere v • n < 0. Finally, the shift operator Sq , when acting on 
v, rotates the velocity through the angle — kip, where ip is the scattering angle (from v' to v). 

When < t < To, [t/To] = 0, and / G (x, v, t) — /i(x, v, t). No recollisions are yet possible and equation (|l|) reduces 
to the standard Boltzmann equation 



D 
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/i(x,v,t)=no / dn(vn)[0(vn)6 n + 0(-vn)] /i(x,v,t) =C B / 1 (x,v,f) , (5) 
Js 1 



where C B is the Boltzmann collision operator. 

When t > To, the distribution /i(x, v,t) splits into two parts. With probability exp(— vTq) the charged particle 
continues to perform free cyclotron motion, having explored the whole circle during the first period Tq. With prob- 
ability [1 — exp(— i/Zo)] the particle suffers collisions and becomes a wandering particle among the hard disks. Since 
the probability of being a wandering particle is less than unity, the distribution for t > To should be renormalized. 
Hence the need for (||). 

The sum in the generalized equation ([!]), from k = 1 to k = [t/To], takes into account all possible recollision events. 
The collision "now" can be the fc'th recollision, with every recollision having the same scattering angle ip, equal to 
the one of the initial collision at t — kTg. Only this initial collision, with incoming velocity Sq v, is described by the 
Boltzmann operator C B , all subsequent ones follow from dynamics alone, weighted by the survival probability from 
one collision to the next, exp(— ^To). Thus, for a sequence of k recollisions on gets the factor exp(— kvTo). Finally, 
summation over k from 1 to [t/To] takes all possible recollision sequences into account. 

On the basis of an intuitive derivation, as sketched above, and by analogy to the results for the general Lorentz 
model, Refs. [^,^| assume that the generalized Boltzmann equation (^) gives an exact description of the time evolution 
of our system in the Grad limit, 

j. _ ( a — > and N/A = n — > oo , . 

Grad \ TIQ = COnSt . 

According to (|^) the radius of the scattering disks approaches zero but their number density increases at the same 
time in such a way that the mean free path of the wandering particle, A = (2na)~ 1 remains constant. 



In the next section, we discuss the initial value problem for the Liouville equation. In Section [II, the effect of 



recollisions is analyzed, and Section IV contains the derivation of the generalized Boltzmann equation. We conclude 



with a short summary. The shift operator along the trajectory involving recollisions is described in the appendix. 



II. DYNAMICAL EQUATIONS 



We denote by T(x, v, t;yi, ...,yjv) the joint probability density for finding the moving point charge at time t at 
point x with velocity v among N hard disk scatterers located at positions yi, y2,...,yjv. F satisfies the normalization 
condition 

/ dv dx dy 1 ...dy N F(-K,v,t;y 1 ,y 2 ,...,yN) = 1 , (7) 
and is the solution of the Liouville equation 
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i ? (x,v,t;y 1 ,y 2 , ...,yjv) = a^Tix- y j ,v)F(x,v,t;y 1 ,y 2 , ...,yjv) 



where T(x — yj, v) is the binary collision operator for scatterer j, 

T(x - y, , v) = / dn (v • n) [0(v ■ n)b n + 0(-v ■ n)] <5(x - y, - an) 
Js 1 



(8) 



(9) 



The operator 6 n acting on the velocity vector has been defined in (^), and an is the point of impact with respect to 
the center of the disk. 

At time t = 0, the initial condition for the Liouville equation (0) is assumed to be of the form 



F(x,v,t = 0;yi, ...,y N ) = /i(x,v,0) 



0(|x- yi |-a) 



.4=1 



p(yi,--,yAf) 

1 — 7rna 2 



(10) 



Here /i(x, v, 0) is the initial state of the particle, and p(yi, ...,yjv) describes the static, non-overlapping distribution 
of scattering disks, 



p( yil ...,y N )= 11 e(\ Yl -y } \-2a)/M , 

l<i<j<N 



(11) 



where Af is the normalizing factor such that J nN dyi ■ ■ ■ yjvp(yi, Yn) = 1- 

To complete the statement of the problem, we shall for simplicity assume that the region Q enclosing the system has 
finite area A and no boundaries (e.g., a two-dimensional torus). We furthermore assume that the linear dimensions 
of O are very large compared to the cyclotron radius. 

The factor [1 — irna 2 ]^ 1 assures the proper normalization (0) of the initial condition since 



/ dyi---ayjvp(yi,---,yjv) 

Jn N 



N 

n 



l* - y<i - a) 



= 1 — nna 



(12) 



Clearly, the probability density /i(x, v, t) is obtained from F(x, v, t; yi, yjv) by integration over the positions of 
all scatterers, 

/i(x,v,t) = / dyi...dy N F(x,v,t;yi,...,y N ) . (13) 
Jn N 

Integrating the Liouville equation (Q) we thus find 

■^/i(x,v,t)=oJv/ dyiT(x-yi,v)/2(x,v,t;yi), (14) 

where 

/2(x,v,<;yi) = / dy 2 ...dy N F{x,v,t]yi,... 1 y N ) . (15) 

Equation ( |l4| ) will be the main object of our study. Using the explicit form of the collision operator (JsJ) , we rewrite 
it as 



D f 

— /i(x,v,t) = an / dn(v-n)[0(v-n)&n + 0(-v-n)J 7 a (x,v,n,t) , 

Ut Jgi 



where 



J a (x, v,n,i) = A / dYF(x,v,t;x - an, Y) . 

J!!"- 1 

where from now on we use the short-hand notation Y = (y 2 , ...yjv)- 



(16) 



(17) 
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In Eq. (|17|), yi = x — an, which corresponds to the precollisional presence of the particle at the surface of the 
scatterer. As has already been mentioned in the Introduction, only velocities such that (v ■ n) < appear as argument 
in the distribution F. 

Our aim here is to demonstrate that Eq. ( |l6| ) becomes a closed kinetic equation for the distribution /i(x, v, t) in 
the Grad limit (||), and that taking this limit we shall recover the generalized Boltzmann equation (Q). 
We note that the Liouville equation (||) implies the relation 

F(x = y + an, v, t; y, Y) = F(S^ (y + an), S^v, t = 0; y, Y) (18) 

where S^P denotes the backwards shift operator along the exact phase trajectory of the iV-scatterer problem. A 
crucial point is to separate clearly two possibilities: (a) The particle did not collide with the scatterer at y before 
time t. (b) It did (leading to the problem of recollisions). The concequences of this distinction will be discussed in 
the next section. 



III. RECOLLISIONS 



For any fixed configuration of scatterers (y, Y) = (y,y2, ••• ) yjv)) one can construct the past history for the moving 
particle leading to the position y + an and velocity v at time t. The probability weight for the initial conditions at 
t = is defined in Eq. @. 

At t > 0, when x = y + an, the particle is about to collide with the scatterer at y, since (v • n) < 0. Depending on 
the past history, we separate the integration domain Q N ~ 1 in ( |i"7| ) into a union of disjoint time-dependent subdomains 
Afc(y + an, v, t) defined as 

Ao(y + an, v, t) = {Y : no collision with scatterer at y before t} , (19) 

and 

A&(y + an, v, t) = {Y : k collisions with scatterer at y before t} , (20) 
where fe = l,2,...,oo. Clearly, we have 

Aj(y + an,v,i) n A fe (y + an, v, i) =0 (21) 

when j ^ k, and 

oo 

|J A k (y + an,v,t) = il"- 1 . (22) 

k=0 

Thus, we may rewrite Eq. ( |l7| ) as 

J Q (y + an,v,n,t) = ^A / dY F(y + an, v, t; y, Y) . (23) 

k=o -' A fc(y+ an . v > t ) 

We now introduce a further partition of each Afc(y + an, v, t), k — 1, .., oo, into two disjoint parts, 

A fe (y + an, v, t) = A<° } (y + an, v, t) U A ^ (y + an, v, t) , (24) 



Af (y + an, v, t) n A« (y + an, v, t) = 0, (25) 

such that 

A^ (y + an, v, t) = [Y : no collisions with scatterers at Y within the time interval [t — kT a , t]} . (26) 

Here T a is the period between two successive collisions of the particle with the scatterer at y. 

If Y 6 A^°' ) (y+an, v, t), then the particle collided with the scatterer at y k times during the time interval [t— kT a , t]. 
Furthermore, it had not collided with this scatterer before t — kT a . We may therefore conclude that 
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YeA[ (,) (y + an,v,t) -F(y + an, v, t; y, Y) = F(S~ k (y + an), S~ k v, t — kT a ; y, Y) (27) 

where S a is the shift operator along the trajectory of the particle between two collisions with the sactterer at y, in 
the absence of other scatterers. We construct this operator explicitly in the Appendix, see also ||. 

In the Grad limit, the volume of Y-space for which collisions with other scatterers occur between recollisions with 
the same scatterer is negligible compared to the volume corresponding to a sequence of consecutive recollisions, and 
in fact 

Urn || AjU ||=0, (28) 

Grad 

where || • • • || represents the relative volume with respect to the set f2 w_1 . The reason behind Eq. ( |2l| ) is that the 
probability to return to the scatterer at y after an intermediate collision with another scatterer involves an extra power 
of a and thus vanishes in the Grad limit. The corresponding problem in the absence of a magnetic field is discussed 
in 0. Close to the Grad limit, (28) allows us to change the volume of integration in (p3|) from A^(y + an, v, i) to 



Any + an,v,t). 

Let us introduce the limit function 

T k (y, % fc v, t - kT ) = Urn ( - ^ ; — ] dY F(y + an, v, t; y, Y) ; A^ = A . (29) 



Grad I l|A^ 0) (y + an,V,i)||y J A W (y+an.v.t) 



where ( p7j ) has been used for convenience. This limit function no longer depends on the number of scatterers, N. 
In terms of Tt, the Grad limit of (K3]) becomes 



It/To] 

7 (y,v,n,t) = l|Ai 0) (y,v,i)|| T k (y, S^ k v, t - kT ) . (30) 

k=0 

Note, with reference to the Appendix, that the shift operator Sq depends on the vector n, even in the Grad limit. 
The relation 

AjjfV.v.t) = Ao(y,S - fe v,i-fcTo)nr(i-fcT ,i) , (31) 

is valid in the Grad limit where the set r(ti, t^) corresponds to those regions of Y-space in which the particle suffers 
no collisions with the N — 1 scatterers at Y during the time interval (t\,t2). 

We now evaluate the volumes of the sets Ao(y, S^ k w, t — kTo) and T(ti,t2). In order to do so, we define the set 
complementary to A (y, v, t), 

Ao(y, v, t) = {Y : at least one collision with scatterer at y before t} , (32) 

If t > To, two types of histories consistent with the condition in j3^ ) are possible: (a) No collisions with other 
scatterers occured within the time interval (i — To,t). (b) The particle collided with other scatterers after the previous 
collision with the scatterer at y. Close to the Grad limit possibility (b) is improbable and in the limit, possibility (a) 
is realized with probability 1. If, on the other hand, t < Tq, the condition in ( |32| ) can only be realized by possibility 
(b), the probability of which vanishes in the Grad limit. Thus, the characteristic function \ is, in this limit, 

_ / ifO<t<T , . 

^A (y,v,t) - | XT(t-T ,t) i£t>T . [66) 

The characteristic function for the set Ao(y, v, t) is 

XA„(y,v,t) = 1 - XA (y,v,t) ■ ( 34 ) 

We now evaluate the characteristic function of the set T(ti, t%). Let C be the trajectory of the particle for t\ < r < i 2 . 
Then, 

£ = {x(r)|*i <r<< 2 }. (35) 
The distance between C and a given point a is 
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dist(£,a)= min llx(r) - all . (36) 

tl<T<t 2 

The characteristic function Xr(t u t 2 ) — f° r finite N — is given by 

N 

Xr(t u t 2 ) = liO[diBt(C, yj )-a}. (37) 
The volume of the set T(t\, t2} is then the integral 

-, N-l 

(38) 



||r(*i,* 2 )|| = 34l^rr / _ dYxnnM- J I dye(dist(£,y)- a ) 



We have on the right-hand side of this expression assumed that we are close to the Grad limit so that the positions 
of the scatterers become independent of each other. We evaluate the integral 

2aL(t 1 ,t 2 )' 



j jf dy0(dist(£,y) - a) 



1 



A 



+ 0{a z ) , (39) 



where L(t%, £2) is the length of the curve C. Combining ( p8[ ) and (p9|), and going to the Grad limit, we find that 

||r(ti,*2)|| = e- 2naL ^^ = e ~ v[ - t2 - tl] , (40) 

where we have used that the particle moves with constant speed v so that L{t\,t2) = vfa — t±). 
Using Eq. @ combined with Eqs. @ and (p3), we find that 



II A I fMI J 1 ifO<t<T 

||A (y,v,t)|| = | 1 _ e _ u2b . fi>To (41) 

We are now in the position to calculate the volume of A^ (y, v, t) through Eq. ( |3l"| ) combined with Eqs. ( f40| ) and 
©• 

The set T(t — kTo,t), consists of all Y such that the particle has k consecutive collisions with the scatterer at y 
without collisions with other scatterers. If we follow the motion of the particle during the time interval 
[t—(k + 1)Tq, t — kT ] that preceeds the time interval for which T(t — kT , t) is defined, and we assume t large enough 
so that t — (k + 1)Tq > 0, two histories are possible: (a) The particle has collided with other scatterers within the 
time interval [t — (k + l)To,t — kTo], or (b) it has collided with the scatterer at y, in which case it has not collided 
with any other scatterer during this period. If (a) is the case, then the configuration belongs to the set A^ ' (y, v, t) 
— see Eq. (|3l| ) — otherwise the configuration belongs to the set T(t — (k + 1)Tq, t). Thus, we have that 

r(t - kT , t) = A^ 0) (y, v, t) U r(t - (k + 1)T , t) . (42) 

Furthermore, 

A( o) (y,v,t)nr(t-(fc+l)T o ,t) = 0. (43) 
The volume of the set A? (y, v, t) is therefore 

||Ai 0) (y,v,i)|| =||r(t-ATo,*)||-||r(t-(fc + l)T ,t)|| . (44) 
Eq. @), which is exact in the Grad limit, gives 

||A[ 0) (y,v,t)|| = e -^ T " (1 - e -" T °) if t>(k + l)T . (45) 

If we now assume that kT Q < t < (k + 1)T so that T(t - (k + 1)T , t) = 0, we have that A ( ° ] (y, v, t) = T(t - kT , t) 
and consequently, 

||Ai 0) (y,v,i)|| =e~ vkT ° if fcT <i<(fc + l)T . (46) 
Using Eqs. ( |45| ) and (|4^) (remember that depends on n), we may rewrite Eq. ( |3C| ) 

I (y, v, n, t) = (1 - e-^o) * £ ' ( y , ^v, t - kT ) + e-»W T °W F [t/To] (y, S- [t/T ° ] v, t 

k=0 ^ 

(47) 

Our next task is to find the relation between /i(x, v, t) and ^fe(y, v, t). 



' t ' 
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IV. THE GENERALIZED BOLTZMANN EQUATION 

The Grad limit of eq. ( |l6| ) becomes 

:/i(x,v,t) =an [ dn(v • n)[0(v ■ n)b n + 0(-v • n)] / (x, v, n, t) , (48) 
> Js 1 



D 

Hi* 



with Iq given by eq. ([47]) with y = x. 

If t < T , as a consequence of Eq. ( p7|) T)( x , v , n , t) = Fo(x,v,t). To express Fo(x, v, i) through /i(x, v, t), we 
consider the right hand side of (|2^) with k = and set for the moment 

F(x, v, t; y, Y) = F W (x, v, t; y, Y) (49) 

to stress that it is a solution of the A^-scatterer problem. Near the Grad limit, i.e., for sufficiently large N and small 
a provided Na — constant, we have 

(x, v, * = 0; y, Y) » I F^ 1 ) (x, v, i = 0; Y) . (50) 

(See Eq. ( ^0| ) for a — > 0.) Assuming that Y G Ao(y + an, v,t) (i.e., no collisions with the scatterer at y before time 
t), we obtain from (pi) 



FW (y + an, v, t; y, Y) = F^ (5^- X) (y + an), S^v, t - 0; yY) (51) 
Combining the last two equalities, we get 

F W (y + an, v, t; y, Y) » I F^" 1 ) (y + an, v, t; Y) . (52) 



for Y e A (y + an, v, t). Thus, near the Grad limit we have from Eq. (|29|) 

F (y,v,^)^ 1 — / dYFC JV - 1 )(y + an ) v ! t;Y) (53) 

||A (y + On,V,t)|| y Ao (y+an,v,t) 

for all t > 0. The simplest case < f < To, leads to the approximate equalities 

A (y,v,<)«^\ (54) 



||Ao(y,v,t)|| «1 



F (y,v,t) w /i(y,v,i) , 



(55) 



which become exact in the Grad limit — see Eq. (|4l|). Hence, for < t < To, Eq. (|4q ) is identical to the usual 
Boltzmann equation (@) for /i(x, v, t). We remark that this is perhaps the simplest way to extract the Boltzmann 
equation directly from the Liouville equations in the Grad limit. 

When t > T , Q N — which is the space over which FW is averaged in order to obtain f\ — is split into two 
disjoint sets, 



(56) 



where 



Aq = {(yi 



,Yn) 



no collisions during the period (t — To, t)} 



(57) 



Both subsets Ao and Ai are defined for any given phase-time point (x, v, t), t > T). If the particle is not scattered 
during the time (t — To, t), there was no earlier scattering due to the periodicity of the motion. Thus, the subsets Ao 
and Ai do not actually depend on time t > Tq and can be equally well defined at t — Tq. The solution of the Liouville 
equation (§]) with initial conditions ([To|) for (yi, Y) G Ao is 



FW(x s v s t;y 1 ,... ) y JV ) = /i(S (-*)x,5o(-t)v,0) 



n »(i5o(-*)x- yi i-o) 



nna^ 



(58) 
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where So (At) is the shift operator along the collisionless trajectory of the particle from t to t + At. Using Eq. J4C|), 
we hnd that in the Grad limit 



HA||-e-" T « , (59) 

and consequently 

\\Ai\\=l-e-^ ■ (60) 
We may write the distribution /i(x, v, t) for t > To as the sum of F averaged over Aq and A\, 

A(x,v,t) = p ||/ 1 ( ° ) (x,v,t) + ||A||/ 1 (1) (x,v,t) , (61) 

where 

/W( XjVji ) = 1 f dydYFW(x,v,t;y,Y) , * = 0, 1 . (62) 
ll-^WII J a, 

We generalize the notation ( |6l|) to < t < To by setting Ao = and A\ = fl N in such a case. 

Let us consider Eqs. fl2^), (p9|), (S) and (^) for t > T . If fc > 1, in (p9|), we may repeat the argumentation leading 
to equality ( p3| ) and obtain 

^(y,5 - fe v,t-fcr )«— ^ ?— / dYf^- 1 )(S -*(y + <m ) J S' -*v,t-Aro;Y) J 

HA^Cy + oSoVv,*)!! A<°)(S- fc (y+a„),v,t) 

(63) 

where A; = 0, [i/Tb], near the Grad limit. A minor change of notation (set N — N + 1 and = y-i-i, i = 2, ...) 
leads to formulas similar to the definition (^2|) of f± with the only difference that the function F is averaged not 

over the whole "ergodic" domain Ax, but over a subdomain A^ C Ai- The final step is to assume that in the Grad 
limit 

My,v,t) =/ 1 (1) (y,v,t) , fc = 0,l... , t>0. (64) 

This assumption means that we can neglect in this limit the difference between average values of T(x, v, t; yi, ...yjv) 
taken over different subsets A^ of the "ergodic" set Ai- This assumption is in fact one of the central postulates 
of kinetic theory (the average behavior of the test particle is the same for almost all configurations of scatterers) 
and we admit it here without proof. The corresponding rigorous result for the linear Boltzmann equation (without 
magnetic field) was obtained in 0. Finally, we note that the representation of /i(x, v, t) in the form ( |6l| ) (cycling 
and wandering particles, see Refs. j|,[|) makes sense only for t > Tq. At the first stage of the motion, < t < To, all 
particles should be considered as wandering, i.e., having a chance to undergo a collision in the future. This explains 
why we put Ai = Q, N and therefore f± = fi for < t < To, sec comment after Eq. ([32|). The above formulas (|4§|), 
( |58| ) — (|62| ) define completely the non-mar kovian kinetic equation (^8| ) in the Grad limit. The final step is to reduce 
the equation to a more convenient form. To this end we define f G by Eq. (0). Thus, we have the normalization 



J dxdv/ G (x,v,t) = | j 



if < t < T 
e- vr ° if t >T . 



(65) 



Eq. (ETj) may thus be written 



[t/T ] 

I (x, v, t) = Y, e^ kTo f G (y, 5 - fc v, t - fcTo) . (66) 

k=Q 



The "non-colliding" distribution function /} ^ satisfies the transport equation 



^/ 1 (o) (x,v,t) = 0. (67) 



The "colliding" distribution function, satisfies the equation 
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— (l- e -"T \ f {1) 

Dt v 



fP(x,v,t)=n / dn(v-n)[0(v-n)6 n + 0(-v-n)] 



[*/t ; 



£ e-^/ G (y,5 - fe v,t-fcT ) 



(68) 



When < > T j we use the definition fcl) and write this equation 



£ / G (x,v,t) = n £ 



[*/To] 



/ dn(vn)[0(vn)& n + 0(-vn)] / G (y, S^v, i - kT Q ) 
Js 1 



(69) 



When t < To, we add Eqs. (|67]) and ( pg) . The resulting equation has the same form as Eq. ( |69| ) — which is the 
generalized Boltzmann equation of Ref. MM . 



V. CONCLUSION 



We have presented a systematic derivation of the generalized Boltzmann equation (Q) from the Liouville equation 
in the Grad limit. The derivation was done at the traditional physicist level of rigor and we stressed several points 
which were not proved mathematically. Moreover, a convergence of the right hand side of Eq. ( |i"4| ) for /i(x, v, t) 
is obviously not equivalent to convergence of the "true" function f[ N ^ to the solution of the generalized Boltzmann 
equation. Therefore, there still remains a considerable effort from a mathematical point of view to prove the validity 
of the generalized Boltzmann equation at the same level of rigor as it was done for the usual Boltzmann equation (see 
Ref. [^| and references therein). However, some of the basic ideas needed and some hints at the difficulties that need 
to be overcome to construct such a proof can be found in the above "physical derivation." 

Another open problem is to clarify the connection of the generalized Boltzmann equation with the BBGKY hierarchy. 
The arguments discussed above show that the validity of the generalized Boltzmann equation follows directly from N- 
scatterer dynamics in the Grad limit. The generalized Boltzmann equation is clearly the non-mar kovian equation for 
a system with strong long-time pair correlations in contradistinction to the usual Boltzmann equation. It is therefore 
important to understand what the correlation functions look like in this case. 

In order to do this, one needs to consider the whole hierarchy and find its asymptotic solution in the Grad limit. 
We hope to do this in a forthcoming paper. 
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VI. APPENDIX 



In this Appendix, we construct the shift operator S a which was defined in (27). It shifts the particle from a position 
right before a collision with the scatterer at y along a collisionless trajectory to a new precollisional position on the 
surface of the scatterer. 

We consider in the following a succession of collisions between the charged particle and a single scatterer. Assume 
that the particle follows the trajectory marked "1" in Fig. [j] This trajectory is a circle arc centered at a distance A 
from the center of the scatterer. The scattering occurs on the surface of the disk, with a scattering angle ip, and the 
particle shifts to a new cyclotron orbit "2". Clearly — (v' • n) = (v • n) = |v| sin(i/>/2), where v' is the precollisional 
velocity and v the postcollisional one. As a result of this symmetry, the center of the new orbit "2" has the same 
distance A from the center of the scatterer as the previous orbit. Orbit "2" leads back to the surface of the scatterer 
where a new collision occurs, and the particle is shifted to a third cyclotron orbit "3" . The direction of n in the first 
collision and n of the second collision is shifted by an angle 2/3, where (3 is given by 

fl A 2 - i? 2 + a 2 

cos/3= — . 70 

2a/\ 

Here R is the cyclotron radius, and a is the radius of the scatterer. 

The angle between two subsequent points of collision is 2/3 as seen from the center of the scatterer. Similarly, the 
angle between these two points, as seen from the center of the cyclotron orbit one to the other, is 27, with 
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R 2 + A 2 -a 2 

C ° S7= 2RA ■ (71) 

The scattering angle in the second collision remains equal to that of the first, if>. Note that in the Grad limit, 
/3 = (it — can take any value from to 7r/2. On the other hand, 7 — > for all collisions in the Grad limit. 

Thus, if we set v = (v,ip v ) and n = (l,ip n ) in polar coordinates, the shift operator S a is, with reference to the 
figure, 

S a (v, (p v ) = (v, f v + ip ~ 27) G£ 5 d (v, (fv+ip) , (72) 

with 7 given by (|^). Similarly, 

S a (l,<p n ) = (1,^ + 2/3) , (73) 

where j3 is given by (p70|) - 

In polar coordinates, the generalized Boltzmann equation (tih becomes 



[t/To] r7I . 

2 



[/ G (x,^, - (k + l)t/j,t-kT ) - / G (x,^, - fe^f-fcTo)] . 

(74) 



This was the form of the generalized Boltzmann equation used in Rcfs. 
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FIG. 1. The charged particle follows orbit "1" and collides with the hard disk, thus shifting to orbit "2". It then recollides 
with the disk and shifts to orbit "3" . The cyclotron radius is R and the distance between the centers of the cyclotron orbits to 
the center of the scatterer is A. The angles separating two subsequent collisions is 2/3. 
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